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We perform near ultraviolet ACS prism spectroscopy of 11 Small Magellanic 
Cloud (SMC) planetary nebulae (PNe) with the main aim of deriving the abun- 
dance of carbon. The analysis of the ACS spectra provide reliable atomic carbon 
abundances for all but a couple of our targets; ionic C 2+ abundances are calcu- 
lated for all target PNe. With the present paper we more than double the number 
of SMC PNe with known carbon abundances, providing a good database to study 
the elemental evolution in low- and intermediate-mass stars at low metallicity. 
We study carbon abundances of Magellanic Cloud PNe in the framework of stel- 
lar evolution models and the elemental yields. Constraining SMC and LMC 
stellar evolutionary models is now possible with the present data, through the 
comparison of the final yields calculated and the CNO abundances observed. We 
found that SMC PNe are almost exclusively carbon rich, and that for the most 
part they have not undergone the hot-bottom burning phase, contrary to ~half 
of the studied LMC PNe. The yields from stellar evolutionary models with LMC 
and SMC metallicities broadly agree with the observations. In particular, evo- 
lutionary yields for M to < 3.5 M Q well encompass the abundances of round and 
elliptical PNe in the SMC. We found that the carbon emission lines are major 
coolants for SMC PNe, more so than in their LMC counterpar ts, indicating that 



metal licity has an effect on the physics of PNe, as predicted by lStanghellini et al. 



fl2003f ) 



Subject headings: Planetary nebulae, stellar evolution, nucleosynthesis, ultravio- 
let spectroscopy, Magellanic Clouds 



1. Introduction 



Planetary Nebulae (PNe) have been studied for decades, and through those efforts an 
understanding has come of the final phases of stellar evolution of stars with masses in the ~1 
- 8 M range (also called the low- and intermediate-mass stars, or LIMS). Planetary nebulae 
are ejected toward the end of the LIMS evolution, during the final thermal pulses on the 
Asymptotic Giant Branch (TP-AGB), thus they are ideal probes to test the ISM enrichment 
from these stars in a quantitative way. The scientific importance of studying PNe in the 
Magellanic Clouds (LMC, SMC) can be very simply summarized: LMC and SMC PNe are 
absolute probes of stellar and nebular brightness, size, and shape, because their distance 
is relatively well known, and because their interstellar extinction is comparatively low and 
uniform. In particular, the SMC is a unique laboratory to study resolved stellar populations 
in a low metallicity environment: chemical analysis of the PNe therein gives a direct test 
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of cosmic recycling at low metallicity, which is essential to understanding the recycling of 
elements within galaxies at high redshift. 

The gas ejected at the PN phase contains elements that have been produced in the 
stellar interiors, and then carried to the stellar surface by the convective dredge-up processes 
(Iben & Renzini 1983, van den Hoek & Groenewegen 1997). The single-star evolution models 
predict a chemical enrichment of the outer region of evolved LIMS, both for Galactic and 
Magellanic Cloud LIMS alike, which can be summarized as follows: 

During the first Red Giant phase, the convective envelope penetrates regions that are 
partially CNO-processed. This dredge-up results in a 13 C and 14 N enhancement, and a 
decrease of 12 C. Afterward, He-burning starts in the stellar core. Later, the hydrogen and 
helium burning occur alternately in two nuclear burning shells, surrounding the CO core, 
and the star ramps up to the AGB. 

The second convective dredge-up (for M to , the turnoff mass, larger than ~3 M ) occurs 
at the onset of the AGB phase, when the H-burning shell is temporarily extinguished. This 
process carries 4 He, 13 C, and 14 N-rich material to the stellar surface. 

During the TP- AGB phase, the envelope is able to dredge-up material after each thermal 
pulse, carrying 4 He, 12 C, and other relatively light s-process elements to the surface. This 
process is called the third dredge-up. 

For M to >3-5 M (exact mass depends on metallicity, Marigo 2001), an additional 
process is thought to occur that alters the chemical composition: the so-called hot-bottom 
burning (HBB) that processes most of the carbon into nitrogen, occurring during the quies- 
cent interpulse periods between thermal pulses. 

The key to assessing the above predictions is to measure the abundances of the pro- 
cessed elements, particularly C, N, and O. Carbon and nitrogen enrichment depends on the 
progenitor mass, yielding to a direct connection between observed progenitor mass (i. e., 
population) and chemical content. By measuring the C and N abundances in PNe, one can 
at once validate key elements of stellar evolution theory, and measure the contribution of 
LIMS to the enrichment of the ISM. While O and N abundance analysis is straightforward 
to do with PNe, owing to their bright optical emission lines, the carbon analysis requires 
satellite UV spectroscopy. 

In our earlier study of LMC PN abundances (Stanghellini et al. 2005) we have shown 
that PN morphology is a surprisingly useful indicator of the progenitor stellar evolution 
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and population. We used a subsample of LMC PNe images available from HST and 
measured the carbon abundances with STIS spectroscopic analysis. We found that nitrogen 
enhancement is correlated with asymmetry. These results are consistent with the predictions 
of stellar evolution only if the progenitors of asymmetric PNe have on average larger masses 
than the progenitors of symmetric PNe. Our results are the first of the kind for extra-galactic 
PNe, and are thus not biased by the large selection effects that limit the observation of PNe 
in the Galactic disk. 

Owing to the smaller metal abundance in the SMC than in the LMC, it is very important 
to extend the LMC study to the SMC. Carbon abundance determination in SMC PNe would 
probe the metal enrichment rates for low metallicity stelar population. 

On these grounds we obtained near UV spectra of 11 SMC PNe that, together with the 
already published LMC and SMC UV data, provide a database that has a statistical impact 
on the study of carbon in Magellanic Cloud PNe. In the present paper we discuss the data 
acquisition and analysis of the ACS prism spectra (§2), the abundance analysis to derive 
carbon in SMC PNe (§3), and a discussion of our results in the framework on PN evolution 
and populations (§4). The conclusions are given in §5. 



2. Observations and Data Analysis 
2.1. Observations 

We observed 11 SMC PNe using the Advanced Camera for Surveys (ACS) prisms 
PR130L in the Solar Blind Channel (SBC) and PR200L in the High Resolution Channel 
(HRC). Our aim is to detect the C IV A 1550, the C III] A 1909, and the C II] A2326 nebular 
emission lines, to determine the ionic abundances of the carbon ions that are excited in the 
PN regime. 

We select our targets from the SMC PNe that have been previously observed by us with 
the Space Telescope Imaging Spectrograph (STIS, Stanghellini et al. 2003; Shaw et al. 2006). 
With this selection we guar antee that their m o rphol o gy, size, opt i cal ex tinction, and optical 



fluxes are known. From the lStanghellini et al.l (120031 ); IShaw et al.l (120061 ) sample we chose all 



PNe that we deemed bright enough to be observed with the observing configuration within 
a few orbits (total H/5 fluxes larger than 2.5xl0~ 14 [erg cm~ 2 s -1 ]). In addition, our selected 
targets have angular sizes smaller than < 0.5" in order to prevent blending of emission lines 



1 LMC and SMC PNe are typically less than 0.5 arcsec across, they can be physically resolved from space. 
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in the slitless spectroscopy. The targets are listed in Table (TJ together with the observing 
log. Of these, a few have been observed before with the IUE, but none has a sound carbon 
abundance determination. 

The SBC detector is a 1024 x 1024 solar-blind Csl Multi-Anode MicroChannel Array 
(MAMA), with ~ 0.034 x 0.030" pixels, and a nominal 35 x 31"field of view. The 
HRC detector is a 1024 x 1024 SITe CCD with ~ 0.028 x 0.025"pixels, covering a nominal 
29 x 26"field of view. The wavelength range of SBC/PR130L is ~ 1200- 2000 A, and the 
useful wavelength range of HRC/PR200L is ~ 1800- 4000 A. The wavelength scale of the 
prisms is non-linear, with spectral resolution decreasing towards longer wavelengths. For 
PR130L, the dispersion varies from about 2 A/pixel at the blue end (R ~ 300), to ~ 10 
A/pixel at 1600 A (R ~ 80) and 30 A/pixel at the red end (R ~ 30 at 2000A). For PR200L, 
the dispersion varies from about 6 A/pixel at the blue en d (R ~ 150), to ~ 20 A /pixel a t 



2500 A (R ~ 60) and ~ 200 A/pixel at 4000 A (R ~ 10) (ILarsen et al.ll2006l : iLarsenl 120061 ) . 



For each target, the HRC observations consist of a direct image through a broad band 
filter (F330W), and two exposures through the prism (PR200L) in order to perform rejection 
of cosmic rays. Similarly, the SBC observations consist of a direct image through F165LP, 
and one exposure of PR130L, unless the exposure time is longer than one orbit, in which case 
two exposures were taken. The direct images have been acquired to establish the zero-point 
of the wavelength scale. 



2.2. Data Analysis 



The pipeline calibration of our data provided bias-su btracted, dark-correcte d, and flat- 
corrected spectral images, as processed with CALACS (jPavlovsky et al.l 120041 ) . We used 
MultiDrizzle and aXe to combine the images, correct for geometric distortion, and to identify 
any bad pixels in the spectrum. The bad pixels were excluded from the analysis of our flux- 
calibrated one-dimensional spectra. 

The position, size, and magnitude of the extracting sources were identified and deter- 
mined using the direct im ages SBC/F165LP and HRC/F330W via analysis with SExtractor 
( iBertin fc Arnoutslll996l ). The position of the extracting sources on the direct images have 
then been projected on the prism images for spectral extraction. 

The spectral extraction was done with the aXe software (iKummel et al. I l2005h in PyRAF@. 
For prism slitless spectroscopy, the trace and wavelength solutions include spatial variations 



2 PyRAF is a product of the Space Telescope Science Institute, which is operated by AURA for NASA. 
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across the HRC and SBC detectors. The wavelength and flux calibrations are provided by 
the ST-ECF group in the configuration files for the aXe software. Two white dwarf stan- 
dards were used to determine the flux calibration for both HRC and SBC prisms. The 
HRC/PR200L wavelength calibration has been secured w ith observations of the LMC PN 
SMP-79 and a quasar with both STIS and the ACS prism ( Larsen et al. 2006h . T he calibra- 



tion of the SBC/PR130L prism was determined using two quasars ( 



Larsen 



20061 ). 



Background subtraction was performed during spectral extraction. To remove the sky 
background, a local background is estimated by interpolating between the adjacent pixels 
on either side of the target spectrum, outside the extracting area. The optimal weighting 
algorithm was chosen to enhance the signal-to-noise ratio of the extracted spectra. The 
algorithm assigns lower weights to pixels which contain only a small fraction of the target 
flux, due to the larger distance from the spectral trace. 

Faint objects near the targets were identified by position and size, and then their con- 
tribution to the flux was subtracted in the extracted spectra. The contaminating flux from 
all other sources was estimated with a gaussian emission model by using the sizes and mag- 
nitudes derived from the direct image. This contamination flux was then subtracted when 
extracting the target spectrum. 

The tabular wavelength and flux of the final extracted spectrum was then read using 
task tprint and converted to 1-d spectrum using the task rspectext. We averaged the 
multiple exposures using task scombine. The final reduced spectra of the 11 PNe are shown 
in Figure [1] to [TTJ The nebular line fluxes have been measured by integrating the area above 
the continuum using the IRAF task splot with the d or w options. 

In Table 2 we give the complete spectral analysis of our targets, including the identifi- 
cation wavelength (column 2), the line ID (column 3), and the mea sured flux as a fraction o f 



Fha (column 4). where the H/3 fluxes are from the HST analysis of IStanghellini et al.l (120031 ) 



and IShaw et al.l (120061 ). It is worth noting that the emission around A1650 could well be a 
blend of H II A1640 and O III] A1663 (the latter corresponding to the 1661/66 A doublet). 
For those low- excitation PNe whose optical spectra do not show the A4686 emission line 
(SMP 18, SMP 20, SMP 24) it is safe to assume that the emission is due to O III]. In the 
other cases we identify the emission as a blend. We do not calculate abundances of other 
elements than carbon in this paper, so we do not present here a detailed analysis of the other 
emission lines. 

Typical uncertainties in the emission line fluxes are ~ 5% in the bright lines. We have 
examined data quality flags of each exposure, and noticed that a few of the PR200L frames 
have some saturated pixels at the red end of their spectra. This is caused by red pile-up, due 
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to the low dispersion in the red wavelengths. In general the saturated pixels do not affect 
the measured fluxes. Only in two PNe, SMP 15 and SMP 18, the saturated pixels affect the 
C II] A2326 emission line, therefore these are listed as lower limits. 

Furthermore, in some cases the spectra are affected by bad pixels. For SMP 24, the 2D 
PR130L spectrum falls into an area of bad pixels, as indicated in the data quality file. If 
we exclude the bad pixels when extracting this spectrum we end up with a very marred ID 
spectrum, thus the C IV A 1550 line measurement is not reliable for this nebula, even if there 
seems to be a feature at the appropriate wavelength. The PR200L spectrum looks much 
better, but the extension shows a dip on the C III] A 1909 line affecting the measured flux 
at an estimated level of 15% level. For SMP 16 the bad pixels in the PR200L frame affect 
the measured fluxes of the C II] A2326 and C III] A1909 lines at an estimated level of 25%. 
Finally, in the case of SMP 18 the C II] A2326 flux might be affected at the 15% level. In 
the case of SMP 25 a prominent NIV] A1485 is seen in the spectrum, but we do not give its 
flux in our tables since there are few columns of bad pixels 1425-1525A range for this PN. 

Ultraviolet prism spectroscopy performed with the ACS have not been acquired often 
in the past, and the analysis modes used in this paper are unique. In order to have a sanity 
check of our calibration and spectral extraction we compare the ACS prism HRC/PR200L 
spectrum of SMP 79, a PN in the LMC, with the corresponding spectrum observed with STIS 
spectroscopy. The prism spectrum was extracted in the same way as for our SMC PN targets. 
We measured the flux of C III] A1909 to compare with the flux measured from the STIS 
spectrum. In the HST archive there are a total of 24 prism exposures of SMP 79 taken at 12 
different positions of the CCD. We extracted the spectrum from each exposure, and measured 
the C III] flux in each spectrum. The average flux is (9.3 ± 0.4) x 10~ 13 [erg cm~ 2 s _1 ], 
compared to the flux 9.4 x 10~ 13 [erg cm" 2 s _1 ] measured from the STIS spectrum. A 
contaminating source happened to lie just at the spectral trace of C II] A2326, so we did not 
attempt to compare this emission line with the one from STIS for LMC SMP 79. Errors in 
our flux measurements are ~ 5% except the specific cases listed above, as marked in Table 
2. 



The flux measurements need to be corrected both for Galactic foreground extinction 
and for the SMC extinction proper. The relation between observed and de-reddened fluxes, 
scaled to H/3, can be written as 



2.3. 



Extinction Correction 




(1) 
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where c is the target- dependent logarithmic extinction at H/3 and fx is the reddening function 
at wavelength A. Since the G alaxy and the SMC have different extinction curves in the UV 
( jHutchings fc Giassonll200ll ). we need to correct separately for each contribution. We can 
write: 

cfx = ccfx,G + csmc/a,smc, (2) 

where the suffixes "G" and "SMC" refer to the Galactic foreground and intrinsic SMC 
extinction, respectively. 

To evaluate the Galactic for eground ex t inctio n we used the Galactic H I foreground col- 
umn density map const ructed by iBqt et al.l (120041 ) from data of the the Parkes H I Survey of 



the Magellanic System (jBriins et al. 
the emission from —60 to +50 km s _1 



20051 ). The column density was obtained by integrating 



outside the SMC velocities. The G alactic color excess 
= 5.8 xfO 21 [atoms cm -2 mag -1 ] ( IBohlin et al.lll978 



map is calculated using N(H I) / Eb-v 
Savage fc Mathislll979l ). Figure [T21 shows the Galactic foreground color excess map, super- 
imposed to the Digital Sky Surve)lj image of SMC. The target-specific foreground Galactic 
extinction constant, cq = 1.47-Eb-V; was estimated from Figure 13 for each target. 

The constant cq is then used to correct the fluxes for foregro und Galactic extinction 
using the extinction function f\G provided by lCardelli et al.l (119891 ): 



— 10 CG ^ A ' G 



(3) 



In order to correct for the SMC extinciton we estimate the optical extinction constant 



as: 



csmc = 2.875 log 



{Ha/H(3\ 
2.85 



(4) 



where Ha/Hf3o is the flux ratio corrected for Galactic foreground. We then correct the flux 
ratios for SMC extinction from: 



h 
h 



— — ) l() l '- ; v< 'a 



l P/ 



(5) 



where the SMC the extinction function is from lPrevot et al.l (J1984J). In Table 2 we list the 
Galactic (column 5) and SMC (column 6) extinciton constants, as calculated above, and the 
final corrected intensity ratios (column 7). Note that the line extinction estimated in the 



3 The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Govern- 
ment grant NAG W-2166. The images of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of these institutions. 
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foreground overestimates the total line extinction in a few cases, thus the SMC extinction 
has been assumed to be zero. 



2.4. Comparison to Previous Observations 



While none of the observed targets had observations in the literature that allowed carbon 
abundance determination, several have UV spectra taken in comparable wavelength ranges, 
and limited comparison could be made among data sets. Let us examine these cases. 



SMP 6: The FOS spectrum published by IVassiliadis et al.l (119961 ) is deemed to be 
inadequate for precise flux measurements (see also SMP 28 below), thus we do not compare 
their fluxes to ours, although we do observe the same bright emission lines. 

SMP 15: A low S/N IU E spectrum of this SMC PN, also known as N 43, has been 
acquired by lAller et al.l ( 19871 ). The only emission line available in the 1200-1910 A wave- 
length window is CIII] A1909, with Fa = 7 x 10~ 13 [erg cm -2 s _1 ], which is within ~ 10 % 
of our measurement. 



SMP 20: Similarly to SMP 15, the only emission line measured by lAller et al.l (119871 ) 
on the IUE spectrum of this SMC PN, also known as N 54, is C III] A1909, with F A = 
7 x 10 -13 [erg cm _2 s _1 ], which is also within ~10% of our measurement. 

SMP 28: The UV spectrum o f this S MC PN short-ward of t he C III] A1909 line is 
available from Meatheringham et al. ( 1990 ). Meatheringham et al. ( 199ol ) give flux values, 
and quote the flux error to better than 15% for brightest lines. By comparing the flux of 
the bright emission lines that we have in common, O IV] A1404, N IV] A 1487, He II A1640, 
and N III A1755, we found that they agree to better than the quoted errors for the two bluer 
lines, but the agreement gets worse for the redder lines measured in our PR130L spectrum . 
In addition to the IUE spectrum there is a FOS spectrum of this PN (IVassiliadis et al.lll996l ) 
whose fluxes are very uncertain, deemed to be systematically off by the authors. For this 
reason we will n ot compare our r e sults with the FOS fluxes. We can use the FOS spectrum 
as presented by IVassiliadis et al.l (119961 ) to confirm that they also observe the C III] A1909 
emission line, contrary to what is derived from the IUE spectrum. Furthermore, the emission 
that we obse rve at 1553 A is mor e likely to be associated with C IV than with Ne V, as 
suggested by iMeatheringham et al.l (119901 ) . In summary, the statement that this PN is very 
carbon poor might not be correct in the light of our observations. 



Finally, lldiart et al.l (120071 ) used the IUE final archived spect ra to measu r e the A 1909 
line for several of the PNe in our sample. The A1909 intensities by lldiart et al.l (120071 ) carry 
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much larger uncertainties than ours, due to the low signal-to-noise ratio of the IUE spectra. 
All intensities agree with our data within the uncertainties, except for SMP 6, SMP 13, and 
SMP 18 where the IUE spectra are noisy and the errors quoted by lldiart et all (120071 ) are 
very large. The IUE spectra would not provide the CIV] nor the [C II] fluxes, thus their use 
to abundance analysis can only establish lower limits to the atomic abundances. 



The comparison of the SMP 28 PR130L spectrum with that of iMeatheringham et al. 



(119901 ) gives us confidence that the C IV] line flux that we measured is accurate. Overall 
we do not have errors larger than 5-10% in this part of the spectrum. Furth ermore, the 
compa rison of the C III] emission lines of our PR200L observations with those by lAUer et al. 
(119871 ) confirms that we have reliable measurements for all the carbon lines. 



3. Abundance Analysis 



The ioni c abundances of C + /H + and C 2+ /H + wer e computed using the nebular package 
in STSDAS (jShaw fc Dufourlll995l ; IShaw et al.lll998l ). We used the line intensities of Table 
2, and T e and N e derived from the diagnostic lines measured by Shaw et al. (in preparation. 



hereaf ter SEA), if available, otherwise we used the plasma diagnostics by lLeisy fc Dennefeld 
(120061 ). hereafter LD06. In Table 3 we present the plasma diagnostics used in this paper; 
columns (2), (3), and (4) give respectively the low- and high-excitation T e and the N e 
used, column (5) lists the references for the spectral lines used to calculate the diagnostics, 
and column (6) gives the excitation class (EC) of the nebulae, which have been estimated 
accordingly to Morgan (1984). 

The C 3+ /H + abundances must be derive d from recombination lines, and we did so on 



the basis of approximate relations from lAUerl (119841 ) (Eq. 5.40). The ionic abundances for 
carbon are listed in Table 4, where columns (2), (3), and (4) give respectively the C + , C 2+ 
and C 3+ abundances in terms of hydrogen. 

It is worth noting that, at our spectral resolution, the A2326 emission line could have 
a [O III] component from the A2321 emission. We calculate the volume emissivities for 
C + and 2+ , the latter from the optical emission lines, for all PNe where A2326 has been 
observed. We found that the ratio of the emission volumes between the carbon and oxygen 
transitions is typically ~50 for the electron densities and temperatures of interest, thus the 
[O III] contribution does not affect our C + ionic abundances. 

In order to calculate the total carbon abundances we follow the discussion by Kingsburgh 
& Barlow (1994), also used by LD06. For low excitation PNe where the C + , C 2+ , and C 3+ 
abundances are available, the total carbon abundance can be calculate simply by summing 
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the ionic contribution. This is the case of SMP 13, SMP 15, SMP 16, SMP 18, and SMP 20, 
whose spectra do not show He II emission thus correction for C +4 would be unnecessary. 
This is also the case of SMP 6, whose emission around 1640-1663 A should be due almost 
entirely to the [O III] component given the PN low optical excitation. 

Next we need to calculate the ionization correction factors (ICFs) for the unseen ion- 
ization stages. We correct for the unseen C + lines in SMP 8, SMP 24, SMP 25, SMP 26 
and SMP 28 with Eqs. All and A13 (KB94). The ionic oxygen abundances used to obtain 
these ICFs have been calculated using the oxygen emission lines in SEA (SMP 8) and LD06, 
the plasma diagnostics of Table 3, and the nebular routines. As it turns out, no correction 
is needed for SMP 8. In SMP 25, SMP 26, and SMP 28 we do not see the C+ emission 
probably due to their high excitation. We still calculate the correction and implement the 
correction, as in LD96, but this has a marginal effect on the final carbon abundances for these 
PNe. We also calculated the ICFs to account for the unseen C 4+ emission in the medium- to 
high-excitation PNe. For planetary nebulae SMP 25, SMP 26, and SMP 28, showing He 2+ 
emission and not N 4+ lines, we use Eq. A20 (KB94). The ionic helium abundances needed in 
these cases were calculated with the line intensities from LD06, the prescription of Benjamin 
et al. (1999), and the plasma diagnostics of Table 3. 

The final ICFs are listed in column (5) of table 4, and the corrected abundances ex- 
pressed in terms of A(C)=log(C/H) + 12 are given in column (6), together with a conservative 
estimate of their uncertainties. 

The larger factor of uncertainty for the total carbon abundances is the uncertainty in the 
electron temperature. The errors listed in Table 4 include the propagation of the uncertainty 
in the T e for all ions, if the plasma diagnostics was available in SEA. For the other nebulae 
we estimate similar errors in the T c and propagate the errors, as given in Table 4, unless 
there are other, larger uncertainties, as discussed below. 

As reported in Table 2 the C 2+ fluxes of SMP 15 and SMP 18 are lower limits; we 
estimate that the total carbon abundances calculated for these two PNe might be underes- 
timated by ~ 5%. 

Finally, while most of the emission lines intensities have errors ~5%, they go up to 15 
and 25% respectively in SMP 24 and SMP 16. These errors are propagated through the 
carbon abundance analysis. 

By considering the sample of this paper and the other carbon abundances in the litera- 
ture, including only those who consider all the present excitation levels (Aller 1987; LD06), 
we found that < C/H >smc=(3.75±3.64) xlO -4 . This is ~1.5 times higher than the same 
average for the LMC (< C/H > LM c=(2.49±2.18) x 10~ 4 , Stanghellini et al. 2005). 
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4. Magellanic Cloud PNe and the stellar evolution models 



Carbon, the fourth most abundant element in the Universe (IClaytonll2003l ). is vigorously 
produced in LIMS, thus it probes their evolution. The a-elements (oxygen, neon, argon and 
sulfur), on the other hand, are produced by nucleosynthesis of Type II supernova and provide 
information about the original composition of the PN progenitor at the time of birth. Oxygen 
may be brought up to the LIMS surface by the third dredge-up, and its abundance should 
be used cautiously in this capacity. 

We use the SMC PN carbon abundances determined in this paper, and those available 
in the literature, and relate them to the abundance of the a-elements, to assess the models 
of stellar evolution and the theoretical yields. In order to compare populations of different 
metallicity we a lso include in the plots the results from our study of carbon abundances in 
the LMC PNe Jstanghellini et al.l l2005h . The abundances of N, O, and Ne for LMC PNe 



used here are from lAUer et al.l (119871 ). and Leis y fc Dennefeld (199 6 , 2006 ). Corresponding 
values of N, O, and Ne in SMC PNe come from lLeisy fc Dennefeld! (119961 ). and SEA. 



In Fig. [13] we show the carbon vs. oxygen abundance^ of PNe in the SMC (filled 
symbols) and the LMC (open symbols). Their morphologies are indicated by symbols of 
various shapes. Small symbols are used for PNe of unknown morphologies, i. e., not yet 
observed with HST or, in a couple of cases, where the morphological class was too uncertain 
to be assigned. Abundances are in the usual scale of A(X)=log(X/H)+12. The PNe separate 
into two groups: The first group are those with C/0<1 (below the line), whose shapes are 
almost exclusively bipolar or bipolar core (BC0), as already noted in Stanghellini et al. 
(2007); these PNe are associated with relatively young stellar populations, such as the disk 
of the Milky Way. And the second group are those with C/0>1, corresponding to the less 
massive progenitors. Most of the SMC PNe lie in the C/0>1 part of this plot. 

The predominance of C-rich PNe in the SMC is consistent with the theoretical expec- 
tation that the third dredge- up is predicted to be favored at lower metallicity, i. e., to occur 
in stars of lower masses and with a higher efficiency (e. g., Karakas, Lattanzio, & Pols 2002; 
Stancliffe et al. 2005). Carbon is converted to nitrogen by hot-bottom burning (HBB). HBB 
is likely active only in relatively high-mass AGB stars (Mgg B ~ 3.5-5 M , depending on 



4 For the sake of clarity we do not plot carbon abundance error bars, that are given in Table 4. Note that 
uncertainties of ^10% would be within the symbol size. 

5 PNe with BC morphologies are those with pinched or barrel-shaped cores similar to known bipolars 
with protruding lobes. This classification is somewhat ambiguous and may include some PNe with unusual 
shapes, and some misclassified elliptical PNe. 
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metallicity ) with very deep convection and hot cores (Marigo 2001). Star formation has 
been far more active in the gas-rich LMC than in the gas-depleted SMC. Possibly, many of 
the PNe that have formed in the LMC have higher average masses than those in the SMC, 
affecting the fraction of high-end mass LIMS that will go through the HBB phase. That 
is, the segregation of PNe with bipolar and BC symmetries below the C=0 line is fully 
consistent with model calculations. It is worth recalling, in fact, that Villaver (2004) found 
a lack of intermediate-mass central stars in the SMC, present in the LMC. 

SMC SMP 22 and SMC SMP 25 are two exceptions to this scenario. They are the two 
morphologically unclassified SMC PNe in the lower-left of the plot: low metallicity (LD06) 
and C/O <1. Not much is known about SMP 22. SMP 25 is quite unique among SMC 
PNe. It is the only PN in the SMC where oxygen- rich dust has been detected (Stanghellini 
et al. 2007), while all other SMC PNe observed with Spitzer/IRS have carbon-rich dust. 
Furthermore, the central star of SMP 25 is much more massive (M~0.82 M Q , Villaver et 
al. 2004) than the typical central in Magellanic Cloud PNe (Villaver et al. 2007), and it is 
located in the eastern region of the SMC facing the LMC, a region that contai ns younger 
and more metal rich clusters (iCrowl et al.l l200ll ) and where iDopita et al.l (119851 ) found the 
kinematically younger PNe to be concentrated. SMP 25 may have undergone recent HBB 
activity. 

In Figure [TU where we show the relation between C/O and the neon abundance, we 
confirm the carbon-abundance segregation of SMC PNe. There is no evidence that neon's 
abundance has been modified in any way, nor that the C/O ratio is modified differently in 
the SMC than in the LMC depending on the initial neon. 

The log(0/N) vs. log(C/N) plot (Fig. [15]) also shows some obvious and interesting 
groupings of data. The segregation of bipolar nebulae on the left side of the plotted line is 
expected of course, since figures 13 and 15 are not entirely independent. However, it is clear 
that the bipolar that exhibit O > C also have N>C: that is, they are "N rich". Symbiotic 
stars and novae stars (not shown; Nussbaumer et al. 1988) are also found in the same region. 
On the other hand, all of the sample with C/O > 1 from Fig. 13 are also nitrogen poor; that 
is, N < O and N< C. Most of the SMC PNe fall into the latter group, of course. Carbon 
stars (not shown; Nussbaumer et al. 1988) are also found in this region. 

Carbon stars probably have not undergone any HBB processing. We can compare the 
data in Fig. 15 with the yields f orm models of stellar evolution. We limit our comparison 
to the fi nal yields calculated by iMarigol (120011) . t o avoid overcrowding of th e plot. The 



yields by Ivan den Hoek fc Groenewegen 



19971 ) or iKarakas fc Lattanziol (120071 ) would look 



very similar to the ones by Marigo! ( 2001 ) on this plot. We indicate the yields from LIMS 



evolution with starred symbols: the four-pointed stars correspond Z=0.008 while the six- 
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pointed stars are for Z=0.004, and we used smaller symbols for stars with progenitor masses 
too low to have undergone the HBB (Mg^ n B - 3.5 and 4.0 M Q for the SMC, and the LMC 
models respectively, as in Marigo 2001). The abundance ratios measured for round and 
elliptical PNe correspond very well with the final yields predicted for the evolution of the 
lower mass-end of LIMS in both the SMC and the LMC. Similarly, the massive models are 
in good agreement with the abundance ratios of C, N, and O in the PNe in the LMC with 
bipolar and BC morphologies. These results are extremely satisfying verifications of difficult 
model predictions. In Figure [TBI we plot the sum of carbon, oxygen, and nitrogen abundances 
vs. carbon of LMC and SMC PNe. Figfl6]is a classical diagram for assessing the efficiency of 
CNO-cycling . The CNO abundance would be preserved if these elements were only acting 
as catalysts, while carbon enhancement via third dredge up is clearly present in the SMC 
PNe. 



Based on a study of LMC and SMC optical lines, IStanghellini et al.l (120031 ) found that 



the observed intensity ratio I[oTn]A5on7/lHfl seem s to increase with metallicity. The pho 



toionization models in IStanghellini et al.l (120031 ) show that metallicity impacts the relative 



emission line strength of the major coolants, and that the UV carbon emission lines have 
a significant effect in cooling PNe with low metallicities. Now that we have observations 
of carbon emission lines in both SMC and LMC PNe we can assess the validity of those 
models. In Fig{T7|we plot the histogram of the IA5007/IA1909 line intensity ratio for the SMC 
(thick line) and the LMC (thin, shaded histogram) PNe. The distributions, which have been 
normalized for the number of PNe in each sample, appear to be different at very low and very 
high intensity ratios, indicating that the carbon emission line is a very important coolant for 
most SMC PNe. We calculate that the average line intensity ratio is 6.23 and 2.36 respec- 
tively for the LMC and SMC PNe, and the median values of Ia5007/Iai909 are 3.34±2.98 and 
1.83±0.83 for the LMC and the SMC PNe respectively, marking a sharp difference in the 
dominant cooling agents at different metallicities. The carbon A 1909 line is a much more 
efficient coolant, with respect to the A5007 line, in the SMC than the LMC PNe, The CIV 
line can also be an important coolant in the low-metallicity, highly ionized PNe. 



5. Conclusions 

The present paper provides a database of 11 SMC PNe whose carbon abundance has 
been accurately derived through prism spectroscopy with the ACS / HST. The number of 
reliable SMC PNe with well-measured carbon abundances has doubled, finally allowing sta- 
tistically useful analysis of carbon production and CNO-cycling in PNe at the very low 
metallicity in the progenitors of SMC PNe. 
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By comparing the abundances of SMC and LMC PNe we found that most observed SMC 
PNe are carbon rich, except in a couple of unusual cases, SMP 22 and SMP 25, where there 
is indication that the progenitors underwent the HBB process. This seems to indicate that 
most SMC PNe derive from low-mass (M to <3.5 M Q ) and low-metallicity progenitors. By 
comparison, the LMC has a much varied PN population, where both carbon-rich and carbon- 
depleted PNe are present, indicating a larger range of PN progenitor mass and metallicity, 
including several PNe whose progenitors could be in the 3-5 M mass range. While a larger 
SMC PN sample would improve the impact of these findings, it is worth recalling that both 
LMC and SMC samples were selected homogeneously, and that with the present observations 
we have increased the statistical significance of the SMC PN sample to almost the same level 
of confidence as the LMC sample (Shaw et al. 2006). 

The CNO abundances in Magellanic Cloud PNe can now be used to test the predictions 
of models of stellar evolution. We found that the data agree impressively with the final 
yields from Marigo (2001) models, calculated both for the SMC and LMC metallicities. In 
particular, the yields calculated for turnoff mass <3.5, and 4.0, M seem to reproduce PN 
abundances of round and elliptical PNe in the SMC, and in the LMC, respectively. Yields 
from the more massive LIMS encompass well the bipolar PNe abundances. If we accept the 
model results then the production of carbon through the CNO cycle in the SMC PNe shows 
that most SMC PN progenitors are in the low-mass end of the LIMS mass range, while the 
LMC PNe could have progenitors in the whole LIMS range. In some cases extra-mixing or 
other process that were not included in the current models are needed to justify the observed 
CNO abundances. 

Compared to the SMC, the LMC has a been a site of recent or ongoing star formation 
and heavy-element enrichment. One of the most important results of our study is that the 
average abundance of carbon in SMC PNe is ~ 1.5 times higher than the average carbon 
measured in LMC PNe. Using the Magellanic Cloud sample of PNe with determined carbon 
abundances, and whose IRS/Spitzer spectra allows a determination of the dust chemistry 
(Stanghellini et al. 2007) we find that <C/H>=(4.69 ±3.33) xl0~ 4 for Magellanic Cloud 
PNe with carbon-rich dust. This is ~35 times greater than the same ratio for PN with 
oxygen-rich dust, <C/H>=(1.35 ±0.88) xlO -5 , confirming a strong correlation between 
dust and gas chemistry. 

The present observations also show that the high-excitation carbon emission lines are 
major or dominant coolants in the low metallicity SMC PNe, as predicted by Stanghellini et 
al. (2003), much more so than in LMC PNe, where the C III] A1909 line intensity is typically 
a small fraction of the [O III] A5007 line strength. 
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Fig. 1. — The extracted spectra from the PR130L (left panel) and PR200L (right panel) 
prism observations in SMP 6. 




Fig. 2.- 



As in Figure 1, for SMP 8. 
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Fig. 3.— As in Figure 1, for SMP 13. 




Fig. 4.- 



As in Figure 1, for SMP 15. 
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6.— As in Figure 1, for SMP 18. 
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8.— As in Figure 1, for SMP 24. 
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Fig. 10.— As in Figure 1, for SMP 26. 
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Fig. 11.— As in Figure 1, for SMP 28. 
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Right Ascension (J2000) 



Fig. 12. — The contours represent the Galactic foreground color excess map towards the 
SMC are overlapped on the the Digital Sky Survey image of SMC. The contours represent 
Eb-v between 0.04 m and 0.08 m from top to bottom, plotted for every 0.01 m interval, with 
the 0.04 and 0.08 contours marked in thick lines. The positions of our targets are indicated 
with crosses. 
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Fig. 13. — The distribution of SMC (filled symbols) and LMC (open symbols) PNe on the 
log(C/H) - log (O/H) plane. Symbol shapes represent the different morphologies. Circles: 
Round, diamonds: Elliptical, triangles: Bipolar core, squares: Bipolar, pentagons: Point- 
symmetric PNe. Small circles refer to unknown morphology, see text. The thick solid line 
represents the 1:1 correlation. The solar symbol indicates the abundances from Asplund et 
al. (2005). 
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Fig. 14.— The distribution of SMC and LMC PNe on the log(C/0) - log (Ne/H) plane. 
Symbols are used as in Fig. Q31 
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Fig. 15.— The distribution of SMC and LMC PNe on the log(C/N) - log(0/N) plane. 
Symbols for data a re used a s in F ig. [T3l The starred symbols represent the yields from 
stellar evolution by iMarigol (120011 ). where the four-points stars are for the LMC and six- 
points stars are for the SMC PNe. Within the same set of models, different points indicate 
different initial mass, and the smaller symbols represent the yields from the evolution of 
M to <Mgg B . The sun locus would be at (0.61; 0.88) in this plot, but it is not marked to 
avoid overlap with a target and few model points. 
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Fig. 16. — SMC and LMC PN abundance sum of C, N, and vs. carbon abundance. 
Symbols are as in Fig. [131 
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Fig. 17. — Cooling of SMC and LMC PNe through the major emission lines. The Ia5007/Iai909 
histogram is shown for the SMC (heavy line) and the LMC (light shaded histogram) PNe. 
Note the very different distributions for Ia5007/Iai909 < 1 an d >6. 
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Table 1 — Continued 
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Table 2. Spectral analysis 
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O^Ofi 

zozo 


P TTl 
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Q1 Q7 C 

yo.y / 


n nq 
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n nn 
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1 no so 

lUZ.oZ 


QA/TP 1 8 
olvlr lo 


100U 


P TAT' 

U IV 
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U.Oo 


n n7fi 
U.U /0 


n nc;s 
U.UOo 


1.1 / 




1 fifi^ 


n TTT 
v^J 111 


Qfi 
z.yo 


n n7R 

U.U / u 


U.UOo 


4 an 
4.yu 




1 QDQ 

iyuy 


P TTTl 
Kj 111J 


1 1 q sn 
no.oU 


n n7fi 

U.U/O 


n n^s 

U.UOo 


1 7Q on 
1 It). ZD 




O^Ofi 

zozo 


P TTl 
Kj 11J 


■~>fiO 1 fi 
-^X>Z. 10 


n n7fi 

U.U/O 


n n^s 

U.UOo 


~^QO Q7 


SMP 20 


1550 


C IV b 


6.88 


0.00 


0.00 


6.88 




1663 


OIII] 


5.93 


0.00 


0.00 


5.93 




1909 


era] 


528.30 


0.00 


0.00 


528.3 




2326 


CII] 


87.06 


0.00 


0.00 


87.06 


SMP 24 


1663 


OIII] 


2.18 


0.047 


0.00 


2.49 




1909 


era] 


147.60 d 


0.047 


0.00 


168.43 


SMP 25 


1397-1407 


[O IV] 


65.93 


0.069 


0.042 


109.10 




1550 


CIV 


68.02 


0.069 


0.042 


104.40 




1640-1663 


He II ,0 III] 


337.30 


0.069 


0.042 


504.40 




1750 


N IIll 


106.50 


0.069 


0.042 


155.10 
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Table 2 — Continued 



Name A ID F A a c G c SM c lA a 



SMP 26 



SMP 28 



1640- 



1909 


cm] 


31.06 





069 





042 


44.87 


1550 


CIV 


509.50 





075 





190 


1853.00 


1663 


He II, III] 


158.20 





075 





190 


524.10 


1909 


cm] 


703.40 





075 





190 


1895.00 


1407 


[0 IV] 


36.48 





098 





041 


65.49 


1487 


[N IV] 


162.00 





098 





041 


275.40 


1550 


CIV 


61.75 





098 





041 


102.20 


1663 


He II ,0 III 


257.30 





098 





041 


413.50 


1750 


N III] 


88.24 





098 





041 


138.00 


1909 


cm] 


50.10 





098 





041 


78.34 



a In terms of H/3=100, where optical fluxes are from lStanghellini et al.l (120031 ). IShawl (120061 ). 

b P-Cygni absorbed component was also detected corresponding to this feature. The flux 
corresponds to the redshifted emission line. 

c Flux uncertainty~ 25% 

d Flux unceratinty~ 15% 
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Table 3. Plasma Diagnostics 



Name T c a T c b N e b ref EC 

[K] [K] [cm-*] 



SMP 


6 


15600 


14400 


11400 


LD06 


4 


SMP 
SMP 


8 

13 




13700+}}^ 
12800™ 


2770 
2900 


SEA 
SEA 


2-3 

4 


SMP 


15 


16200 


12000 


5000: 


LD06 


2-4 


SMP 


16 




11800 


5000: 


LD06 





SMP 


18 




11860+*°° 


3590 


SEA 


0.5 


SMP 
SMP 


20 
24 




13820:™ 
11620^ 


3880 
2780 


SEA 
SEA 


1-2 
1-2 


SMP 


25 




21100 


9800 


LD06 


6-7 


SMP 


26 




18000 


300 


LD06 


8 


SMP 


28 




19200 


20700 


LD06 


8 



a Temperature for low excitation ions 
b Adopted by Authors in reference 
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Table 4. Carbon Abundances 



Name 


C+/H+ 




C 2 +/H+ 




C 3+ /H+ 




ICF(C) 


A(C) 


SMP 6 


1.016x10" 


-5 


1.917x10- 


-4 


2.184x10- 


-5 


1.0 


8-35±8:S 


SMP 8 






9.596x10- 


-5 


3.631x10- 


-5 


1.0 


o 1 9+0.17 

^-J-^^o.is 


SMP 13 


2.877x10- 


-5 


4.753x10- 


-4 


3.445x10- 


-5 


1.0 


o 7Q+0.05 
°- '-J-0.15 


SMP 15 


> 4.889 x 10" 


-6 


1.621x10- 


-4 


1.379x10" 


-5 


1.0 


8.26 + °;°7 


SMP 16 


7.357x10- 


-5 


7.761x10- 


-5 


4.403x10- 


-6 


1.0 


8.19 + °; 08 


SMP 18 


> 6.422 x 10- 


-5 


1.377x10- 


-4 


7.047x10- 


-7 


1.0 


o o-i +0.05 
°- oi -0.20 


SMP 20 


2.694x10- 


-5 


1.492x10" 


-4 


1.283x10- 


-6 


1.0 


8 2 c i +0 - 2 


SMP 24 






1.485x10" 


-4 






1.02 


o -i o+0.05 
S-J-S-0.25 


SMP 25 






1.626x10- 


-6 


1.602x10- 


-6 


1.38 


6 64 +0 - 09 

u -°^-0.13 


SMP 26 






1.349x10- 


-4 


6.497x10- 


-5 


1.45 


o Afi+0.04 
<5- 4b -0.05 


SMP 28 






4.255x10- 


-6 


2.525x10- 


-6 


1.34 


6 Q6+ - 02 
°- yu -0.03 



